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DETECTION  CAPABILITY  OF  LINEAR-AND-POWER  PROCESSOR 
FOR  RANDOM  BURST  SIGNALS  OF  UNKNOWN  LOCATION 

INTRODUCTION 

There  is  always  a  need  for  reliable  detection  of  weak  signals 
in  noise.  However,  the  problem  is  greatly  confounded  when 
knowledge  of  the  signal  characteristics  or  structure  is  minimal 
or  nonexistent.  Yet,  this  is  frequently  the  situation  in  initial 
detection  of  a  distant  echo  from  an  unknown  target(s),  when 
virtually  nothing  is  known  a  priori. 

An  example  is  furnished  by  attempting  to  detect  weak  targets 
at  long  range  by  means  of  high-frequency  underwater  acoustic 
transmissions.  In  addition  to  the  deleterious  effects  of  the 
background  noise  and  the  reverberation  (due  to  surface/bottom 
reflections  and/or  suspended  sediments  for  example),  the  signal 
echo  (if  present)  is  greatly  attenuated  in  strength  due  to  the 
distance  of  propagation.  Furthermore,  the  echo  signal  generally 
possesses  a  random  structure  due  to  the  nature  of  the  target, 
medium,  and  boundaries. 

Classification  and  identification  of  target  echoes  cannot  be 
reliably  accomplished  without  first  achieving  reliable  detection. 
For  weak  echoes  especially,  it.  is  mandatory  to  know  that  the 
signal  processing  techniques  being  employed  are  operating  near 
the  absolute  limits  possible  in  the  particular  environments 
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encountered.  For  random  signals  with  little  structure  or 
location  information,  these  limits  were  not  known  and  have  just 
begun  to  be  established  (see  references  1  through  6).  This 
report  continues  this  investigation  for  burst-like  signals. 

It  is  therefore  necessary  to  formulate  and  solve  the  problem 
of  optimum  or  near-optimum  detection  of  random  signals  with 
little  structure.  In  particular,  for  signals  that  have  bursts  of 
energy  in  time,  but  of  unknown  locations  and/or  durations  as  well 
as  strength,  the  likelihood  ratio  processor,  which  maximizes  the 
detection  probability  for  a  specified  (tolerable)  false  alarm 
probability  P^,  must  be  derived.  Alternatively,  for  signals  that 
occupy  globular  regions  in  time-frequency  space  of  unknown  shape, 
size,  and/or  location,  the  likelihood  ratio  and  the  attendant 
signal  processing  indicated  by  the  likelihood  ratio  test 
statistic  must  be  derived. 

Since  the  optimum  likelihood  ratio  processor  frequently 
cannot  be  realized  practically  —  because  it  requires  knowledge  of 
unknown  signal  parameters,  one  must  be  satisfied  with  a 
suboptimum  processor,  hopefully  derived  as  a  physically 
reasonable  approximation  to  the  optimum  processor.  However, 
acquiescing  to  a  suboptimum  processor  necessarily  degrades  the 
performance  that  can  be  attained.  Therefore,  it  is  imperative  to 
evaluate  the  extent  of  degradation  that  accompanies  the 
particular  suboptimum  processor  adopted  and  ascertain  its 
acceptability. 
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Several  major  technical  problems  must  be  addressed  and 
solved.  First  is  the  actual  derivation  of  the  likelihood  ratio 
processor  under  conditions  that  are  realistic,  including 
specifically  the  burst-like  structure  of  the  echo  signal  strength 
in  time  and/or  frequency.  The  next  task  will  require 
manipulating  the  exact  likelihood  ratio  test  into  an  approximate 
form  or  forms  that  eliminate  any  dependence  on  signal  parameters 
that  will  not  be  known  in  practical  applications.  For  example, 
the  signal  strength  might  be  presumed  known  during  the  derivation 
of  the  optimum  likelihood  ratio  processor,  only  to  be  eliminated 
in  determining  a  practical  suboptimum  processor. 

The  exact  level  of  performance  attainable  by  the  suboptimum 
processor  is  of  paramount  importance  and  must  be  evaluated. 
However,  since  it  is  suboptimum,  there  will  definitely  be  a  loss 
of  performance  relative  to  the  likelihood  ratio  processor,  which 
loss  must  be  quantified.  If  this  loss  in  performance  cannot  be 
ascertained  exactly  by  analysis,  it  can  be  determined  by  using  a 
combination  of  simulation  and  bounding  procedures.  Direct 
comparisons  of  the  receiver  operating  characteristics  (ROCs)  are 
the  best  means  by  which  useful  quantitative  measures  of  loss  can 
be  furnished. 

The  particular  detection  problem  that  will  be  investigated 
here  is  that  of  a  contiguous  burst  signal  of  known  duration  or 
extent,  but  of  unknown  strength  and  location,  within  a  specified 
search  region  of  given  size. 
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OPTIMUM  PROCESSOR  FOR  UNKNOWN  SIGNAL  LOCATION 

The  problem  will  be  couched  in  the  time  domain  but  could  be 
extended  easily  to  frequency  or  combined  time-frequency  space. 

A  burst  of  energy  is  emitted  and  reflected  off  a  distant 
distributed  target  of  interest.  The  received  weak  signal  echo 
(if  present)  is  accompanied  by  wideband  noise  that  occupies  the 
entire  time  scale.  Since  the  target  distance  is  unknown,  the 
received  signal  time  location  is  unknown. 

It  is  presumed  that  the  received  noise  has  been  normalized  to 
unit  level  (through  past  observations)  and  that  a  search  region 
of  N  time  bins  has  been  established.  The  width  of  each  time  bin 
is  related  to  the  inverse  bandwidth  of  the  received  signal  burst, 
while  the  search  size  N  depends  on  the  uncertainty  of  the  target 
range.  It  is  presumed  that  the  received  waveform  is  passed 
through  a  passband  filter  centered  on  the  transmitted  carrier 
frequency,  the  output  of  which  is  squared-envelope  detected  and 
sampled  sequentially  in  time.  For  an  unknown  Doppler  shift,  it 
will  be  necessary  to  broaden  the  filter  passband. 

DERIVATION  OF  LIKELIHOOD  RATIO 

The  joint  probability  density  function  (PDF)  of  the  N 
noise-only  bin  outputs,  hypothesis  Hg,  is 

Pg (Ui , . . . ,uN)  =  exp(-u^)  •••  exp(-uN)  for  each  un  >  0  .  (1) 
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For  signal  present,  hypothesis  ,  it  will  be  assumed  that  the 
signal  duration,  that  is,  the  number  of  occupied  bins,  consists 
of  M  contiguous  time  slots,  where  M  is  known .  This  corresponds 
to  a  continuously  distributed  target  in  range  with  random 
properties.  For  equal  average  signal  powers  S  in  each  of  the  M 
occupied  bins,  the  conditional  joint  PDF  of  the  N  observations, 
when  the  initial  occupied  signal  bin  is  number  m,  is 

M-l+m _ fa  exp(-au  J') 

crti(ul,  .  .  .  ,uN)  =  pQ(u1,  .  .  .  ,uN)  M  exp(-u  )  '  (2) 

n=m  ^  n 

where  a  s  1/(1  +  S).  Therefore,  the  joint  PDF  of  the  observation 
under  is,  for  equally  likely  initial  starting  points 
m  c  { 1 ,N+1-M} , 

N+l-M  1  M-l+m 

p1(u1, . . . ,uN)  =  p0(Ul, . . . ,uN)  C  k+1^M  TT  *a  exp(wun)}  ,  (3) 

m=l  n=m 


where  weighting  w  s  S/(l  +  S).  Therefore,  the  likelihood  ratio 

for  random  observation  or  measurement  x, is 

1  N 


i(xi 

Vxi 


. ..,x„)  M  N+l-M 

JN  _  3  r - 1 

•••'V  "  N+1*M  5=i 


exp 


r  M-l+m  ' 

W  )  I  X 
L - 1  n 

l  n=m  ) 


(4) 


The  corresponding  likelihood  ratio  test  is  to  compare  the  m-sum 
above  with  a  fixed  threshold,  namely,  in  expanded  explicit  form, 


exp 


+ 


+ 


4* 


exp 


H1 


‘N+l-M 


+  x. 


NJ 


(5) 
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This  same  result  for  the  likelihood  ratio  test  was  derived 
previously  in  reference  4,  pages  19-21.  However,  no  numerical 
results  on  the  performance  of  processors  (4)  or  (5)  were 
presented  at  that  time.  A  connection  to  the  generalized 
likelihood  ratio  test  was  also  made  in  reference  4,  page  22. 

Several  comments  about  optimum  processor  (5)  should  be  made. 
First,  it  is  easily  realized  and  computed  (when  duration  M  and 
weight  w  are  known)  because  it  requires  only  the  evaluation  of 
N+l-M  exponentials,  a  very  reasonable  burden;  this  is  in  direct 
contrast  to  the  case  of  arbitrary  structureless  signals,  where 
(N-M) !  operations  can  be  required.  However,  since  weight 
w  =  S/(l  +  S)  depends  on  the  generally  unknown  signal  level  S, 
processor  (5)  is  not  a  practical  processor.  In  addition,  signal 
extent  M  is  presumed  known,  and  that  knowledge  is  used  in  optimum 
processor  form  (5).  Notice,  also,  how  the  known  contiguous 
signal  structure  is  taken  into  account  by  performing  linear  sums 
of  the  available  data  {xnJ ,  using  the  known  duration  M ,  before 
resorting  to  nonlinear  combinations  over  all  the  equally  likely 
possibilities . 

The  exponential  function  grows  rapidly  with  increasing 
argument.  This  means  that  the  largest  linear  sum  in  equation  (5) 
will  be  significantly  accentuated  and  will  tend  to  dominate  the 
other  exponentials.  This  observation  will  lead  to  a  suboptimum 
processor  with  potential  for  good  signal  detectability. 
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SPECIAL  CASES  OF  OPTIMUM  PROCESSOR 

When  M  =  1  in  equation  (5),  the  likelihood  ratio  test  becomes 

N  V 

YZ  exp(wx  )  v  .  (6) 

n=l 

This  processor  has  already  been  investigated  in  reference  4, 
figures  2  and  12.  Also,  the  signal  has  no  burst  duration  in  this 
case,  since  only  one  signal  bin  is  occupied;  that  is,  there  is  no 
useful  burst  information  when  M  =  1. 


When  M  =  N  in  equation  (5),  the  likelihood  ratio  test  becomes 


N 


n 

n=l 


x 


n 


> 

< 


v  , 


(7) 


where  the  single  exponential  and  the  scaling  w  have  been 
eliminated,  because  they  constitute  a  monotonic  transformation, 
and  do  not  change  the  ROCs.  Processor  (7)  has  been  analyzed  in 
reference  1,  pages  21-22  and  81-90.  Also,  the  signal  location 
has  no  uncertainty  since  it  occupies  the  entire  search  region; 
that  is,  there  is  no  useful  burst  information  when  M  =  N. 


Due  to  the  restricted  behavior  of  these  two  special  cases, 
M  =  1  and  M  =  N,  and  since  they  have  already  been  considered 
previously,  attention  is  confined  here  to  signals  with  a 
nontrivial  burst  structure  with  some  uncertainty  of  location, 
namely,  1  <  M  <  N. 
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APPROXIMATE  LIKELIHOOD  RATIO  PROCESSORS 


If  one  denotes  the  inner  sum  on  n  in  equation  (4)  as  random 
variable  XM(m),  the  likelihood  ratio  test  in  equations  (4)  and 
(5)  can  be  expanded  as 


N+l-M  . 

HZ  exP  (w  XM(m) 

m=l 


N+l-M  00  fj 


m=l  /j=  0 


00  fj  N+l-M 

C  ^  C  <(»)  $  v 

fj= 0  m=l 


(8) 


Thus,  the  fundamental  data  processing  on  observation  {xn}  has 
the  following  typical  component: 


N+l-M 

HZ  xM(m> 

m=l 


N+l-M  (M-l+m 

IZ  HZ 

m=l  l  n=m 


n 


(9) 


due  to  its  form,  processor  (9)  is  called  the  LAP  (linear  and 
power)  processor.  The  attractive  feature  of  form  (9)  is  its 
independence  of  the  generally  unknown  weight  w.  A  reasonable 
approximation  to  the  likelihood  ratio  test  is  furnished  by 
comparison  of  the  statistic  in  equation  (9)  with  a  threshold;  the 
choice  of  control  parameter  n  is  not  obvious  at  this  point, 
although  u  =  1,  2,  3,  4  are  good  candidates  and  should  be 
investigated.  In  fact,  fj  need  not  be  limited  to  integers  in  LAP 
processor  (9).  This  processor  was  previously  suggested  in 
reference  4,  page  21. 


A  limiting  case  of  processor  (9)  is  afforded  by  taking  the 


1/u  root  of  the  left-hand  side,  and  letting  u  +  ».  (The  root 
operation  is  a  monotonic  transformation  and  does  not  change  the 
ROCs.)  The  end  result  is  the  maximum  processor, 


9 


lim 

yU-»°° 


fN+l-M 


_  m=l 


1  /» 


max{XM(m) }  for  m  e  {1,N+1-M}  . 
m 


(10) 


This  maximum  processor  can  also  be  seen  to  be  an  approximation  to 
the  optimum  processor  (4)  or  (5),  since  the  maximum  XM(m)  term  is 
exponentiated  and  tends  to  dominate  the  sum  on  m. 


For  fj  =  1,  equation  (9)  can  be  simplified  to  the  form 
xx  +  2  x2  +•••+  M  xM  +  ••■  +  M  xN+1_M  +  •••  +  2  xN_1  +  xN  >  v  (11) 

(for  M  <  N/2 ) ,  which  indicates  that  the  interior  samples  of  the 
observation  {xnJ  should  be  weighted  more  heavily.  A  similar 
behavior  holds  for  the  quadratic  and  cubic  terms  in  equation  (9). 
This  suggests  the  following  possibility  as  an  alternative 
approximation  to  the  likelihood  ratio  test: 


N 


E 

n=l 


> 

< 


v  . 


(12) 


The  weights  {rn}  should  be  symmetric  and  peak  at 
n  =  (N  +  l)/2;  their  exact  dependence  is  not  obvious,  but  the 
trapezoidal  weighting  in  equation  (11)  is  an  initial  candidate. 
The  choice  of  power  \>  could  be  taken  in  the  range  (1,3)  at  first, 
but  need  not  be  integer.  If  sequence  {rnJ  is  chosen  to  be 
independent  of  M,  equation  (12)  is  a  practical  processor  that 
could  be  applied  in  signal  detection  situations  where  neither  the 
signal  duration  nor  location  are  exactly  known. 
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SPECIAL  CASES  OF  LAP  PROCESSOR  (9) 


When  M  =  1  in  LAP  processor  (9),  the  approximate  test  takes 
the  form 
N 

'  v  .  (13) 

This  test  is  recognized  as  the  standard  power-law  processor, 
which  has  been  thoroughly  analyzed  in  references  3-5,  when 
power-law  \>  there  is  identified  with  power  u  here.  Therefore, 
there  is  no  need  to  investigate  the  M  =  1  case  here  for  different 
values  of  power  /j  in  LAP  processor  (9). 

For  M  =  N  in  LAP  processor  (9),  the  approximate  test  takes 
the  form 

(C  *„)  <  V!  or  II]  xn  <  V  .  (14) 

n=l  n= l 

(The  removal  of  power  n  does  not  alter  the  ROCs.)  Test  (14)  is 
recognized  as  the  standard  energy  detector,  which  has  been 
analyzed  in  references  3-5;  see  the  \>  =  1  case.  Thus,  there  is 
no  need  to  investigate  the  M  =  N  case  here  for  different  values 
of  power  fj  in  LAP  processor  (9).  Furthermore,  the  ROCs  for  the 
M  =  N  test  in  equation  (14)  are  independent  of  u;  therefore,  all 
the  burst-signal  LAP  processors  in  equation  (9)  have  identical 
performances  when  M  =  N. 


E 

m=  1 
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PERFORMANCE  OF  APPROXIMATIONS  AND  OPTIMUM  PROCESSORS 

Processor  forms  (11)  and  (12)  can  be  analyzed  exactly  for  any 
set  of  weights.  Thus,  ROCs  can  be  constructed  fairly  easily  for 
these  two  techniques.  However,  suboptimum  LAP  processor  (9)  will 
require  simulation,  due  to  the  nonlinear  (power-law)  interactions 
that  create  statistically  dependent  random  variables  prior  to  the 
sum  on  m. 

Optimum  processor  (5)  will  definitely  require  simulation  to 
determine  the  ROCs.  However,  it  can  be  simulated  readily  within 
reasonable  execution  time,  since  it  requires  only  N+l-M 
exponentials  per  trial.  This  absolute  upper  bound  on  performance 
requires  and  uses  knowledge  of  M  and  w  =  S/(l  +  S).  It  can  be 
used  to  quantitatively  ascertain  the  losses  that  must  be 
tolerated  when  one  employs  the  various  approximate  processors 
listed  above. 

The  availability  of  an  upper  bound  on  detectability  perform¬ 
ance  allows  one  to  investigate  numerous  weighting  schemes  {rn}  in 
power-law  processor  (12),  in  a  search  for  the  best  possible 
practical  processor;  this  includes  a  search  over  power-law  values 
v.  If  the  losses  between  the  optimum  and  practical  processors 
can  be  reduced  to  an  acceptably  low  level,  by  choices  of  \>  and 
{rn},  the  latter  processors  can  be  employed  with  confidence,  even 
when  there  is  a  lack  of  information  regarding  signal  parameters, 
such  as  location,  duration,  and  strength. 
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OPTIMUM  PROCESSOR  FOR  UNKNOWN  SIGNAL  LOCATION  AND  DURATION 

In  this  section,  in  addition  to  unknown  signal  location,  the 
duration  M  of  the  signal  energy  burst  is  unknown  but  assumed  to 
lie  in  the  range  {A,B} .  In  particular,  M  can  take  on  the  values 
A,  A+l , ,  B  with  probabilities  Q^,  QB»  respectively. 

Furthermore,  for  a  given  duration  M,  the  initial  occupied  signal 
bin  m  can  start  equally  likely  at  one  of  the  positions 
m  =  1,  2,  .  .  .  ,  N+l-M. 

Under  hypothesis  Hq  and  unit  noise  level,  the  joint  PDF  of 
the  observation  is  again  given  by  equation  (1).  Under  hypothesis 
,  if  M  signal  bins  are  occupied,  the  average  signal  level  per 
bin  is  the  same  for  all  the  occupied  bins  and  is  denoted  by  SM< 
This  dependence  of  the  average  signal  level  per  bin  on  M  allows 
for  consideration  of  more  general  cases.  For  example,  the  total 
received  average  signal  level  is  M  if  M  bins  are  occupied; 
this  product  could  be  kept  constant  as  M  varies,  meaning  that 
SM  =  k/M  in  this  particular  case. 

The  joint  PDF  of  the  observation  under  ,  conditioned  on  M 
signal  bins  being  occupied,  is  a  slight  generalization  of 
equation  (3)  to 

N+l-M  .  M-l+m 

P1(u1,...,uN|M)  =  pQ(u1,  .  .  .  ,uN)  XZ  n+1-M  "I  f  *aM  exP(wMUn)}  ' 

m=l  n=m 

(15) 

where  a^  &  1/(1  +  S^) ,  wM  s  SM/(1  +  S^) .  Then,  the  unconditional 
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PDF  is  given  by 


B 


?1 ( U1 '  *  *  ’ ' UN )  _  ^  ^ U1 '  *  *  • ' un I ^ ^ 


M 

B  Q„  a“  N+l— M 


M-l+m 


=  P0(ul' 


x_»  XlV/r  '-*■  XI  I  J.  A  X  -  X  X  I  ALL 

V  C^C  exp  „M  YZ  »n  • 

M=A  m=l  n=m 


(16) 


For  given  data  observation  {xn}  for  1  <  n  <  N,  this  leads  to 
the  likelihood  ratio  in  the  form 


LR  s 


Pl(Xl,...,XN)  Q 


M 

a„  N+l— M 


Pq(Xi  ,  .  .  .  ,xKT) 


N‘ 


M=A  m=l 


(17) 


where  the  linear  sum 


M-l+m 

Vm)  =  YU  x 


n=m 


n 


for  1  <  m  <  N+l— M  ,  A  <  M  <  B  . 


(18) 


Some  observations  on  optimum  processor  (17)  are  in  order  at 
this  point.  This  processor  is  relatively  easily  realized  and 
computed  (when  signal  powers  { }  and  a  priori  probabilities  {QM} 
are  known).  However,  since  these  quantities  are  generally 
unknown,  equation  (17)  is  not  a  practical  processor.  Notice  that 
the  known  contiguous  signal  structure  is  taken  into  account,  by 
performing  linear  sums  (18)  of  the  available  data  {xn}  using  the 
hypothesized  duration  M,  before  resorting  to  weighted,  nonlinear 
combinations  in  equation  (17). 
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To  convert  form  (17)  into  a  practically  useful  processor,  one 
approximates  the  exponential  according  to  exp(y)  =  k  y^,  where 
scale  factor  k  and  power  p  can  be  chosen  to  give  a  reasonable  fit 
over  an  extended  range  of  y.  Then,  equation  (17)  becomes 


B  N+l-M 

LR  =  ^  XM(m) 

M=A  m=l 


(19) 


where 


Q„  a?!  w(J[ 
_  ,  M  M  M 

aM  “  K  N+l-M 


=  k 


M 


5M 


N+l-M 


(1 


+  V 


M +fj 


for  A  <  M  <  B 


(20) 


Implementation  of  processor  (19)  requires  knowledge  of  the 
signal  powers  {SMJ  and  the  prior  probabilities  { QM } *  if  these 
quantities  are  known  approximately  or  estimated,  the  scale 
factors  {aM}  in  equation  (20)  can  be  calculated.  On  the  other 
hand,  if  some  of  this  information  is  absent,  it  may  be  necessary 
to  guess  at  their  behaviors  and  use  a  rough  estimate  for  the 
weights  {o:^}.  If  no  information  is  available,  the  weights  {«MJ 
might  be  taken  as  flat  over  the  range  { A , B } .  In  the  latter  case, 
the  approximate  likelihood  ratio  processor  (19)  becomes 


B  N+l-M 

m  m  <(m)  >  V  .  (21) 

M=A  m=l  x  ' 

The  best  power  /j  to  use  in  tests  (19)  or  (21)  is  yet  to  be 
decided  on;  it  could  well  depend  on  M  and  range  { A , B } ,  as  well  as 
the  average  signal  level,  and  need  not  be  integer. 
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Test  (21)  says  to  perforin  linear  sums  of  the  available  data 
{xn}  over  the  hypothesized  duration  of  length  M,  raise  it  to  a 
power,  sum  over  all  possible  initial  locations  m,  and  then  sum 
over  all  possible  durations  M  within  the  expected  range.  This  is 
a  combination  of  a  linear  processor  and  a  power-law  processor. 

How  well  it  performs  relative  to  optimum  processor  (17)  can  only 
be  determined  through  extensive  simulation.  This  situation  and 
processor  has  not  been  investigated  quantitatively  here;  that 
issue  is  left  for  future  study. 
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DETECTION  PERFORMANCE  OF  PROCESSORS 


The  observed  data  {xn},  1  <  n  <  N,  occupy  a  search  region  of 
size  N  bins,  while  the  signal  (if  present)  occupies  M  contiguous 
bins  of  unknown  location  within  the  search  region.  Three 
processors  will  be  investigated  quantitatively  in  this  section. 
All  require  preliminary  evaluation  of  the  N+l-M  linear  sums 

M-l+m 

XM(m)  =  )  j  x  for  1  <  m  <  N+l-M  . 

M  1 — 1  n  - 

n=m 

The  three  processors  of  interest  are 


Optimum 


(5) 


max{XM(m) }  ^  v  Maximum  (10) 

m 

N+l-M  .  .// 

)  1  [xM(m)J  v  Linear  and  Power  .  (9) 

m=l 

The  power  //  in  the  LAP  processor  is  a  control  parameter  that  can 
be  varied  for  maximum  detectability,  while  S  in  the  optimum 
processor  is  the  average  signal  power  per  occupied  bin.  Because 
the  average  noise  power  per  bin  is  1,  S  is  also  the  signal-to- 
noise  power  ratio  per  bin. 


17 


The  ROCS,  namely,  detection  probability  Pd  versus  false  alarm 
probability  P^,  for  Maximum  processor  (10)  are  presented  (solid 
curves)  in  figures  1  through  9  for  M  =  2,  4,  8,  16,  32,  64,  128, 
256,  512,  respectively.  The  label  S(dB)  on  the  curves  is  equal 
to  10  log(S).  These  simulation  results  were  all  based  on  using 
one  million  independent  trials  of  processor  (10). 

Superposed  on  each  figure  (dashed  line)  is  the  performance  of 
Optimum  processor  (5)  for  one  representative  signal-to-noi se 
ratio;  in  particular,  the  S(dB)  values  investigated  are  9,  7, 
4.25,  2,  0,  -2,  -4,  -6,  -8,  respectively.  These  values  were 
chosen  so  that  the  corresponding  ROC  would  pass  close  to  the 
standard  operating  point  Pf  =  0.001,  Pd  =  0.5.  It  is  immediately 
obvious  from  figures  1  through  9  that  the  corresponding  ROCs  for 
the  Maximum  and  Optimum  processors  are  virtual  overlays  in  every 
case,  there  being  a  very  slight  gain  (<0.1  dB)  in  favor  of  the 
Optimum  processor.  That  is,  the  Maximum  processor  for  detection 
of  a  signal  of  unknown  strength  and  location  is  virtually  optimum 
in  terms  of  its  performance. 

It  was  found  that  the  ROCs  for  LAP  processor  (9)  were 
uniformly  poorer  for  finite  values  of  power  p  than  for  p  =  ®,  the 
Maximum  processor.  Accordingly,  the  ROCs  for  p  =  1,  2,  3,  4  have 
been  relegated  to  the  appendix;  results  for  M  =  2,  4,  8,  16,  32, 
64,  128,  256,  512  are  presented  there. 
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From  figures  1  through  9,  it  is  possible  to  determine  what 
values  of  signal-to-noi se  ratio  S(dB)  are  required  to  exactly 
achieve  the  standard  operating  point,  when  using  the  Maximum 
processor.  When  these  values  are  scaled  by  M,  the  number  of 
occupied  bins,  they  represent  the  total  required  received 
signal-to-noise  ratio.  These  total  signal-to-noise  ratios  are 
plotted  versus  M  in  figure  10  as  the  bottom  curve,  labeled  y  =  ». 
Since  this  Maximum  processor  is  virtually  optimum,  according  to 
figures  1  through  9,  the  bottom  curve  in  figure  10  is  essentially 
an  absolute  lower  bound  on  all  possible  detection  procedures 
functioning  under  the  prescribed  conditions. 

The  remaining  four  solid  curves  in  figure  10  correspond  to 
the  LAP  processor  and  are  labeled  y  -  1,  2,  3,  4.  These  values 
of  required  signal-to-noise  ratio  are  extracted  from  the  ROCs  in 
the  appendix.  They  show  a  steady  degradation  in  performance  as 
power  y  is  decreased  from  y  =  ®;  in  fact,  the  y  -  1  processor, 
which  merely  sums  up  energies  in  a  linear  fashion,  suffers  a  very 
significant  loss  and  should  not  be  used,  except  when  M  is  near  N. 
The  results  in  figure  10  are  somewhat  similar  to  those  for  the 
power-law  processor  in  references  3-5,  where  the  signal  had  no 
structure  whatsoever. 

For  comparison,  superposed  in  figure  10  (dashed  line,  labeled 
OPT-B)  is  the  required  signal-to-noise  ratio  for  the  optimum 
processor  for  detection  of  a  signal  with  no  structure  (see 
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reference  4).  The  difference  in  performance  between  the  OPT-B 
curve  and  the  //  =  »  curve  represents  optimal  usage  of  the  extra 
information  available  in  knowledge  of  the  signal  structure, 
namely,  a  contiguous  burst  in  this  case.  For  M  =  1  and  M  =  N, 
there  is  no  difference  between  these  two  optimum  curves,  because 
there  is  no  extra  burst  information  in  these  cases.  However,  for 
the  intermediate  case  of  M  =  64,  the  additional  burst  information 
leads  to  a  maximum  improvement  of  2.88  dB;  also,  as  can  be 
expected,  the  optimal  processor  form  is  significantly  changed, 
from  a  2.5  power-law  device  (in  the  case  of  no  signal  structure) 
to  a  Maximum  operation  on  partial  linear  sums  (for  burst-like 
signal  structure). 


Figure  1.  ROC  for  Maximum  Processor  with  M  =  2,  N  =  1024 
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Figure  10.  Total  Sigiial-to-Noise  Ratio  Required 
for  Ff  =  0.001,  Pd  -  0.5,  N  =  1024 


25/(26  blank) 


SUMMARY 


Detection  of  a  random  burst  signal  of  known  extent  M,  but  of 
unknown  location  and  signal  strength,  can  be  accomplished  in  a 
near-optimum  fashion  by  means  of  the  Maximum  processor. 
Preliminary  linear  sums  of  duration  M,  starting  at  every  possible 
initial  bin,  are  conducted,  and  the  maximum  sum  is  compared  with 
a  threshold,  for  declaration  of  signal  presence  or  absence. 
Knowledge  of  the  average  received  signal  strength  is  not  required 
or  used.  The  discrepancy  between  the  Maximum  and  Optimum 
processors  is  less  than  0.1  dB  at  the  standard  operating  point 
for  search  size  N  *=  1024  and  for  signal  durations  ranging  from 
M  =  1  to  M  =  1024. 

When  the  signal  extent  M  is  also  unknown,  an  additional 
search  over  the  allowed  range  of  durations  must  be  conducted. 
Simulation  of  the  optimum  (likelihood  ratio)  processor  is 
possible  within  a  reasonable  amount  of  computer  time;  however, 
that  numerical  study  has  not  been  conducted. 
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APPENDIX 

OPERATING  CHARACTERISTICS  FOR  LAP  PROCESSOR 

The  LAP  (linear  and  power)  processor  is  characterized  by  test 

N+l-M  ,  .  y  M-l+m 

5  !  fx^noj  <  v  /  with  =  'EZ  xn  •  ( a— i ) 

m=l  n=m 

The  LAP  receiver  operating  characteristics  for  powers  y  =  1 ,  2, 

3,  4  and  signal  durations  M  =  2,  4,  8,  16,  32,  64,  128,  256,  512 
are  collected  in  this  appendix;  the  search  region  is  fixed  at 
size  N  =  1024  in  all  cases.  The  ROCs  are  based  on  one  million 
independent  trials  of  test  ( A— 1 ) ;  thus,  the  smaller  false  alarm 
probability  estimates  P^  are  not  reliable. 

It  will  be  observed  from  these  characteristics  that  the 
performance  of  LAP  processor  (A-l)  improves  monotonically  with 
power  y ,  regardless  of  the  value  of  M;  in  fact,  the  Maximum 
processor  ( y  -  ®)  is  virtually  optimum,  as  was  shown  in  figures 
1  through  10  of  the  main  text.  The  total  signal-to-noi se  ratios 
in  figure  10  for  y  -  1,  2,  3,  4,  required  to  realize  the  standard 
operating  point  Pf  =  0.001,  Pd  =  0.5,  were  extracted  from  figures 
A-l  through  A-36  here.  Other  higher  quality  operating  points  can 
also  be  investigated  from  the  figures  in  this  appendix. 
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Figure  A-24.  ROC  for  LAP  Processor  with  y 
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Figure  A-35. 


p< 

ROC  for  LAP  Processor 


with  u  =  4,  n  m  256,  N 


p( 

Figure  A-36.  ROC  for  LAP  Processor  with  fj  =  4,  M  =  512,  N 
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